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Abstract. We investigate correlation functions in a periodic box-ball system. For the second 
and the third nearest neighbor correlation functions, we give explicit formulae obtained 
by combinatorial methods. A recursion formula for a specific iV-point functions is also 
presented. 
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1 Introduction 

<3\ \ 

\ A periodic box-ball system (PBBS) is a soliton cellular automaton obtained by ultradiscretizing 

^ ; the KdV equation [TJ [2]. It can also be obtained at the q — >• limit of an integrable lattice 

■ model (a generalised 6 vertex model) Let V\ be a 2-dimensional complex vector space 

O , with basis en and e\. If we consider N tensor product space of Vt, V := Vf N , the transfer 

matrix of the generalised 6 vertex model is a map (endmorphism) T^{x\q) : V —> V with 

a spectral parameter x and a deformation parameter q. The positive integer k denotes that 

the dimension of the auxiliary vector space is k + 1. An important property of the transfer 

^ | matrices is their commutativity: Tfc(u; q)Ti(v; q) = Ti(v; q)Tk(u; q) for arbitrary k, I, u, v. Hence 

they consist a complete set of diagonal operators and the lattice models are integrable. It 

is also noted that the Hamiltonian of the quantum XXZ spin model is essentially given by 

Hxxz '■= nm 7& log^ife <?)• If we take a limit lim T(x; q) =: T, the transfer matrix T maps 
x— >l ° q—>0, x—tl 

each monomial \i1i2 ■ • • ijv) := % e,- L2 ® ■ ■ ■ <g> e*j n € V to a monomial. By identifying \i1i2 ■ ■ ■ *jv) 
with a 10 sequence ii?2 • • • *at, 7" gives the time evolution of the PBBS. Using this relation, we 
can obtain several important properties of the PBBS such as the conserved quantities, a relation 
to the string hypothesis, a fundamental cycle and so on. 

One of the main problems of quantum integrable systems now is to obtain correlation func- 
tions which is fairly difficult even for the XXZ model and the 6 vertex model [6]. Since 
a correlation function of the XXZ model or the 6 vertex model is, roughly speaking, trans- 
formed to a correlation function of the PBBS in the limit q — > 0. For example, a two point 
function (sfsp of the XXZ spin model is transformed to the probability that both ith and 
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|http://www.emis.de/journals/SIGMA/ISG2010.html| 
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Figure 1. Definition of T L for / 6 fix- 



jth componets of the PBBS are 1. Hence, from the view point of integrable lattice models or 
quantum integrable models, it may actually give some new insights into the correlation functions 
of the vertex models themselves to obtain correlation functions of the PBBS. 

In [5], we gave expressions for Appoint functions using the solution for the PBBS expressed in 
terms of the ultradiscrete theta functions. We also gave expressions for the one point, the nearest 
and the second nearest neighbor correlation functions. Note that the correlation functions 
defined there and in this article are not the direct q — > limit of the correlation functions of 
the corresponding lattice model; such a limit generally gives trivial results [5]. We define the 
correlation functions on the phase space of a set of fixed conserved quantities. In this article, 
we show that the second nearest neighbor correlation function can be simplified with some 
combinatorial formulae, and give an explicit formula for the third nearest neighbor correlation 
function by combinatorial methods. A recursion formula for a specific Appoint function is also 
presented. 

The PBBS can be defined in the following way. Let L > 3 and let = { / | / : [L] — > 
{0, 1} such that t}/" 1 ^!}) < L/2 } where [L] = {1, 2, . . . , L}. When / G L is represented as 
a sequence of 0s and Is, we write 



The mapping Tl : Ql — > Ul is defined as follows (see Fig. [JJ : 

1. In the sequence / find a pair of positions n and n+ 1 such that f(n) = 1 and f(n + 1) = 0, 
and mark them; repeat the same procedure until all such pairs are marked. Note that we 
always use the convention that the position n is defined in [L], i.e. n + L = n. 

2. Skipping the marked positions we get a subsequence of /; for this subsequence repeat the 
same process of marking positions, so that we get another marked subsequence. 

3. Repeat part 2 until one obtains a subsequence consisting only of 0s. A typical situation 
is depicted in Fig. [TJ After these preparatory processes, change all values at the marked 
positions simultaneously; One thus obtains the sequence Tj^f. 

The pair (Ql,Tl) is called a PBBS of length L [2j[7]. An element of Ql is called a state, and 
the mapping Tl the time evolution. 

The conserved quantities of the PBBS are defined as follows. Let Qj (/) be the number of 10 
pairs in / marked at jth step in the definition of the mapping Tl. Then we obtain a nonincreasing 
sequence of positive integers, Qj(f) (j = 1,2, ... ,m). This sequence is conserved in time, that 
is, 



For example, (Qi,Q2, Q3, Qa) = (3, 3, 2, 1) for / given in Fig.[TJ As the sequence (Q\, Q 2 , ■ ■ ■ , Q m ) 
is nonincreasing, we can associate a Young diagram with it by regarding Qj as the number of 



/(l)/(2) •••/(£). 



Qj(f) = Qj(T L f)^Qj 



(j = i,2,... 



m). 
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squares in the jth. column of the diagram. The lengths of the rows are also weakly decreasing 
positive integers. Let the distinct row lengths be Pi > P 2 > ■ ■ ■ > P(_ and let rij be the number 
of times that the length Pj appears. The set {Pj, nj}j =l is another expression for the conserved 
quantities of the PBBS. 

First we summarize some useful properties of the PBBS. We say that / has (or that there is) 
a 10-wall at position n if f(n) = 1 and f(n + 1) = 0. Let the number of the 10-walls be s and 
the positions be denoted by a\ < a 2 < • • • < a s . 

We introduce a procedure called 10- insertion. For (61, b 2 , ■ ■ ■ , &l) € Z>o x Z>o x • • • x Z>q (= 

L 

(Z>o) L ) and d := Yl h, the 10-insertion is defined as I(b\, 62, ... , b£) := I 2 (bi, 62, • • • , bi) h, 

i=i 

where I\ : Ql — > &L+2s is the mapping: 

f 1 (n = aj + 2j-l, j G [s]), 



(hf)(n) 



(n = aj + 2j, j € [s]), 

f(n) (l<n<ai), 

/(n - 2j) ( aj + 2j + l<n< a j+1 + 2j, j € [s - 1]), 

f(n-2s) (a s + 2s + 1 < n < L + 2s) 



and the mapping h{bi, 62, ■ ■ • , : 
(J 2 (6i,6 2 ,...,6 i )/ / )(n) 



^L+2(d+s) is defined to be 



f /'(n - 26 J _ 1 ) (n = j + 2b^ h j G [L]), 

1 (n = j + 2Sj-_i + 2/i - 1, j G [L], /i G [6j]), 

(n = i + 2&J--1 + 2/i, j G [L], he [bj]), 



where /' = I\f G Ql+2s, bo = and 6j := G [s] | a« < j} + X^i=i ^i- For example, 
/ = 0011100111000001101000111000000, 

hf = ooiii|io|ooiii|io|oooooii|io|oi|io|oooiii|io|oooooo, 

^(19, 122, 129)/ =OO1110OOlllJ10OOOOO110O10OOOli)1110OOOOlOOO, 

where (lg, 1 2 2 5 129) denotes the 10 sequence whose 9, 22, and 29th elements are 1 and the others 
are 0, that is, 



'J-9; I22, 129) 



(0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0), 



and the expressions |Io| and 10 denote the 10s inserted at / i->- I^f and at I±f h-> / 2 (&i, bi, ■ ■ ., 
bL)(hf) respectively. 



2 Correlation functions of PBBS 

We consider the PBBS with L boxes and M balls. As shown in [5], the correlation functions 
are defined on the states with the same conserved quantities. Let Y be the Young diagram 
corresponding to the partition of M: 

(P 1 ,P 1 ,...,P 1 ,P 2 ,P 2 ,...,P 2 ,...,P e ,P e ,...,P e ), 

ni 712 n l 

t 

where P\ > P 2 > • ■ ■ > Pn > 1 and ^ n jPj = M. The number of squares in the jth largest raw 

i=i 

of Y is Pj. In the PBBS, the conserved quantities are characterized by the Young diagram Y [8]. 
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When we denote by fiy the set of states with conserved quantities given by Y, the iV-point 
function of the PBBS is defined as 



{si,s 2 ,...,s N ) Y = t^— r ^2 /( s 0/( s 2)---/(sjv)- 



\n Y \ 

Since the TV-point function (si, si + di, . . . , s\ + cIn-i}y does not depend on the specific site si 
(because of translational symmetry), we denote 

Cy(d 1 ,d 2 , ■ ■ ■ ,d N -i) = (si,si + di,...,si + djv-i)y, 

where 1 < d\ < d 2 < • • • < djv-i < L. The 1-point function Cy(0) := (si)y and the 2-point 
function Cy(l) := + l)y are easily calculated as Cy(0) = -jj and Cy(l) = ^f^- where 

t 

S = Y, n i- 

i=l 

Although general iV-point functions are difficult to evaluate by elementary combinatorial 
methods, the correlation functions Cy(l, 2, . . . ,N — 1) = (si, si + 1, s\ + 2, . . . , s% + iV — l)y 
enjoy a simple recursion formula. Let us put 

{nj (i = Pj), - ^X, 

fci := > kj, L := L, Lj := - 2fcf, 
U otnerwise, t - : 

for i = 1, 2, . . . , P\. Note that ki and fcj are the number of the rows in Y with the length i and 
that with the length greater than or equal to i respectively. For the Young diagram Y, let Yi be 
the Young diagram corresponding to the partition 



I max [Pi — i, 0] , max [Pi — i, 0] , . . . , max [Pi — i, 0] , 

ni 

max [P 2 — i, 0] , max [P 2 — i, 0] , . . . , max [P 2 — i, 0] , . . . 



Il 2 



, max [P^ — z, 0] , max [P^ — i, 0] , . . . , max [P^ — i, 0] I . 



The following two lemmas immediately follow from the definition of I\ and the procedure of 
10-insertion. 

Lemma 1. In a state f = I\f\ (fi € ^y), no pattern with 101 exists. 

Lemma 2. The number of the pattern 11 in f\ € ily 1 is equal to that of the pattern 111 in I\fi, 
and it does not depend on the specific state f\. If we denote it by M 2 , it is calculated as 

Pi 

N 2 = Y,h{i-2). (1) 

Then we have the following theorem. 
Theorem 1. 

( V1 U...,Hl).(|n, i ;y.Jw...,r) (, = 1,2,...). 
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In particular, 



C '- (1 - 2) = L(W, 1 -l) § t '"- 2) ' 

Proof. We call a 'block' for a repeated 10 pattern of the form 10 10- • -10. We also call a 111 
pattern for the pattern of the form 111. 

First let us consider the case r = 1. From Lemma [21 the number of 111 patterns in I\f\ is 

Pi , L X 

equal to £ h(i - 2). Let B\ be a set {B x £ (Z> ) Ll | E h = h, Bi = (fa, b 2 , . . . , b Ll ) }. In 

i=3 i=l 

the process I\f\ I(Bi)fi, if 10 blocks is inserted between the left 1 and the middle 1 of a 111 
pattern, then the 111 pattern in I(B\)fi disappears. Let : Z>o x (F2 x F2) —> {0, 1} be the 
function: 



^i(b;ui,u 2 ) 



1 (6 = 0, ( Ul ,u 2 ) = (1,1)), 
otherwise 



and tpi(fi) '■= (J{ i € [L\] | (fi(i),fi(i + 1)) = (1,1)}. Since a state of U Y is obtained by 
10-insertion to a state in Q Yl , we have 



Cy (1,2) 



From 



Pi 



L 



L 



L 



Bill fin 



Si 



Si 



Si 



E E E( J (£i)/i)« (j(5 1 )/ 1 )(^+2) 

J BlGBl/ief7y 1 1=1 

Li 

E E E^^aw./i^+i)) 

PieBi /ien Yl i=i 

Li 

E E £*i(&il/i(0,/i(< + l)) 
PieBi/ien^ i=i 

£ 2 ^i(/i)*i(6 i; l,l) 
BieSi/iea ri 

2 ^i(/i) E *i(fc;M) 

/leOyj BieBi 



^Li+fci-2-j / j 



fci 



L|Bi| I |O n | 



E fc(/o 



(L X - 1) 



Pi 



/iGn Yl 



L(Li + fci-l) 



5>(i"2). 



i=3 



Pi 



Hi ~ 2) = ~ 1) " k = L x C Yl (1), 

i=3 i=2 

we find 

Li (Li - 1) 



Cy(l,2) 



L (Li + fci - 1) 
For general r, if we define 

^r(h, fa, . . . ,b r ;ui,u 2 , ■ ■ ■ , u r+ i) 



(2) 
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f 1 ({h,b 2 ,...,b r ) = (0,0,... ,0), (ui,u 2 ,...,u r+ i) = (1,1,... ,1)), 
I otherwise, 

and i'r(fi) ■= U{ i G [£i] | {fi{i),fi(i + !),•• • + r)) = (1, 1, . . . , 1) }, we can proceed in 
a similar manner to the case r = 1 as 

1 L 

Cy(l,2,...,r + 1) = — — — £ £ ••• (/(BOAjCi + r + l) 

1 Ll 
= L | g n^| E E E*-(^-,6*4r;/l(0,-,/l(i + r + l)) 

1 111 ri1 B r eBi /ien n »=i 

E E E^( 6 i'---'^+i'" /i(o,---./i(i+»-+i)) 



1 111 Yl < B r eB! han Yl i=i 

rlgMo I E E Vv(/l)*r(&l,-" ,6r+ll (1,1, ■■■,!)) 
^Li+fci-r-l^ 



X 



r-r E ^(/i) 



LiC yi (l,2,...,r) ^j (Li-i) 



1 ( Ll 

which completes the proof. ■ 

3 The second and the third nearest neighbor 
correlation functions 

3.1 Correlation function CV(2) 

Using 10-insertion and the Young diagram Y\ defined in the previous section, we can rewrite 
C y (2) as 

-H^^E i/ EE«B 1 )/ 1 ) WW B 1 )/ 1 )(n + 2 ). 

From this expression we can prove the following proposition which gives a key formula to cal- 
culate Cy(2). 

Proposition 1. Let <E>o : ^>o x (F2 x F2) -> Z>o be the function defined by 

' (6 = 0), 

b (6/0, ( Ul ,u 2 ) = (0,1), (1,0)), 
k 6-1 (6/0, ( Ul ,u 2 ) = (0,0), (1,1)). 



$ (b;ui,u 2 ) := < 



Then we have 



C Y (2) = j(M 2 + n 2 ), (3) 
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where N% is given in (p} , 
1 

Wi 



BiGBi bG(F 2 ) 2 

and Mh;b) :=tl{*e[Li]|(/i(i),/i(* + l)) =b}. 

Proof. We call a 1 * 1 pattern for the pattern of the form 1*1 where * is either 1 or 0. From 
Lemmas [U and El we already know the number of 1 * 1 patterns in Iif\. Hence we estimate 
the number of the 1*1 patterns generated and eliminated by inserting blocks 10 10 • • • 10 in the 
process I\f\ h4 I(B\)f\. Since there appears no 1 * 1 pattern between blocks, we have only 
to count the number of the patterns 101 appeared in a block and a boundary and that of the 
patterns 111 eliminated by 10-insertion; then, for b n > 1, we observe the number of 1*1 patterns 
appeared by the operation of I2 (B\ ) as 

(1) /i(n)/i(n+l) =00 — > 01010- •• 100 => +(b n - 1), 

(2) A(n)/i(n + l) = 10 — >• 101010- •• 100 =► +b n , 

(3) /i(n)/i(n + l) = 01 — > 010 10- •• 101 +b n , 

(4) / 1 (n)/i(n + l) = ll — »• 11010- ■■ 101 +(b n -l). 

Here we write the referred 1, 0s in I\fi in bold scripts (1, 0) for clarifying the results. Note that 
if fi(n)fi(n + 1) = 11 and b n = 0, then we have the patterns 111 or ll|i~o| or 1110 after the 
operation of ^(-Bi)- Hence, we have 

712 = \B I in 1 ^ ^ t^' 101 ' G /(jBl)/l } " ^' n1 ' G disa PP ears in 
- W ^r l E E i>(fciAM,A(i + i». 

1 111 ri1 BiGBi/ienyj i=i 

Here ${'101' € denotes the number of '101's in I{Bi)f\ and jj{'lll' € I\f\ disappears 

in denotes the number of '101's which disappear by the operation ^(-Bi)- Since 

(61,62, ... ,b Ll ) G Bi — > (6i,6i + i,...,6 il ,6i,...,6i_i) G Si, 

we have 

$ (6 J ;/i(i),/ 1 (i + l))= £ $ (6i;/i«,/i(i + l)). 
B1GB1 B1GB1 

Thus we have 

Li ii 

^ E*o(6i;/i(i),/i(* + l))= E E <I) °( 6 ^W'^ + 1 )) 

= E E fo(/i;&)*o(&i,&). 

BlGBi 6G(F 2 ) 2 

Furthermore, if / and /' are the states of PBBS with the same number of boxes, balls and 
solitons, then c/>o(/;&) = 4>o(f';b) for any b G (F2) 2 > in particular, 0o(/i;b) = (f>o(f[;b) for 
/1 j /{ € ^ Y x • Therefore we obtain ([3]) . ■ 
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For the evaluation of Cy (2), we use the identity for binomial coefficients given in the following 
lemma. 



Lemma 3. For k < I < n, it holds that 



E 

m=0 



k + m 
k 



k + m — 1\ 1 f n — k — m 



k 



Theorem 2. 



C H 2)4(i>( ! -2> + AeLziL + 



I - k 



2k 2 h 



(4) 



L 
1 

I 



U=3 



Li + fci - 2 Li + fci - 1 



Li + fei - 2 Li + fci - 1 



(5) 
(6) 



Proof. There are ( Ll_t fc ^ 1 _ r ^ 1 2 ) elements with b\ = m (m > 1) in the set B\. Since 0o(/i> (0, 1)) 
= o (/i; (1, 0)) = fc 2 and o (/i; (0,0)) + o (/i; (1, 1)) = £i - 2fc 2 , we have 



771=1 



^ ^ o (/i;b)$o(6i;b) = E 

BieBi 6g(F 2 ) 2 
fci 

-E 



7 Li + fei — m — 2 



k\ — m 



Mfi\b)$ (m;b) 

6e(F 2 ) 2 



777=1 



Li + /«i — m — 2 
fci — m 



(m 

6=(0,0), (1,1) 



Li£(m-1)( 

777=1 ^ 



Li + k\ — m — 2 
k\ — m 



+ 2k 2 



b=(0,l), (1,0) 

Li + h - 2 
ki - 1 



that is, 



Cy{2) = - 



fci 

£iX>-i)( 
£^•-2) + —^ — 

7 = 3 



Li + k\ — m — 2 
k\ — m 



2h 



\Bi\ 



+ 



Li + ki-2 
ki - 1 



\Bi\ 



\ 



When we put n = L\ + k\ — 2 and I = L\, turns into 

'Lx + ki-2-2-m 

ym tin 

777=0 ^ 

'L x +ki-l 



Lx+kx-2 
Lx 



ki —2 — m 



) = X>-i; 

' 777=1 



Since \B 



fci 



£(m + l) 

, we obtain ([5]). Equation ([H]) follows from 



L\ + k\ — 2 — m 
fci — m 



3.2 Correlation function Cy(3) 

The evaluation of Cy(3) can be done in analogous way to that of Cy(2), although we have to 
consider various patterns so that the expression may become fairly complicated. 

Using 10-insertion and the Young diagram Y 2 defined in the previous section, we have 



Cy(3) 



1 



L\Bx\\B 2 \\ny 2 \ 



E E E 

SieBi b 2 £B 2 / 2 erv 2 
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Y J {l{.B l )I{B 2 )f 2 ){n){l{B 1 )I{B 2 )f 2 ){n + ^ 



71=1 



where B 2 = { B 2 € (Z> ) 



X] 6i = £; 2 , -B 2 = (61, b 2 , . . . , bi 2 ) } . We call a 1**1 pattern for an 



array of 4 elements of the form 1 * *1 where ** is 00 or 01 or 10 or 11. The following lemmas 
are easily understood by the definition of 10-insertion and Proposition [TJ 

Lemma 4. In a state f = I\f\ (/1 £ £Iy-\)> no pattern with 1101 and 1011 exists. 

Lemma 5. For B 2 € B 2 and f 2 € Qy 2 , number of the 1*1 patterns in fx := I(B 2 )f 2 is 
equal to that of the patterns of the form 1001 and 1111 in I\f\, and is given by 



J^k i (i-3)+j(B 2 



i=A 



where 



j(B 2 ) := ]T $ (& i; / 2 (i), f 2 (i + 1)) (£ 2 = (61,62,---, h 2 ) € B 2 ) 



Lemma 6. Let N 2 := Yl h{i - 3) and TZ 2 : = ^2|(/ 1 ,Bi)->(/ 2 ,B 2 ) = ]ET\ _S jX^)- TVien 



i=4 



have 



B 2 £t3 2 



C Yl (2) = -L(A/- 2 + 7e 2 ). 

^1 



Let us define 



\Bi\\B 2 \\n Y2 \ 



E E E UCllll'/lOOl' ehI(B 2 )f 2 }. 



BieSi B2&B2 / 2 eni 



From Lemmas HHSl we find A/3 = LiCy-^). If we define IZ3 by 

Cy(3) = 2-(A/- 3 + ft 3 ), 

it denotes the term which arises from 10-insertion I(Bi), that is, 
1 



\Bi \\B 2 in 



V2I 



E E E (H'HOl'.'lOH'G/W/^/a} 



S1GS1 B 2 eB 2 f 2 &Q Y2 
- tl{'100r, '1111' € hI{B 2 )f 2 disappear in I(B 1 )I(B 2 )f 2 } 

Then we obtain the following proposition which gives a key formula to calculate Cy(3). 
Proposition 2. Lei us define the functions; $1 : Z>o x (F 2 ) 4 — >• {0, 1,2}: 



$i(6;iti,M 2 ,«3, "4) := < 



'2 (6/0, (m, U2, u 3 , u 4 ) = (0, 1, 1,0), (1,0, 0, 1)) , 

/ («1,U 2 ,«3,«4) = (1,0, 0,0), (0,1, 0,0), 

(1,1, 0,0), (0,0,1,0), (1,0,1,0), (1,1, 1,0), 
' (0,0, 0,1), (0,1, 0,1), (1,1, 0,1), (0,0, 1,1), 
V (1,0, 1,1), (0,1, 1,1) 

otherwise, 



(7) 
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$ 2 :(Z> ) 2 x(F 2 ) 4 ^ {-1,0,1}: 

( 1 (W^O, (ui,U2,U3,«4) = (1,0,0,0),(0,0,0,0)) ; 



$ 2 (.b,b';ui,u 2 ,u 3 ,u 4 ) : = < 

, otherwise 

and $ 3 : (Z> ) 2 x (F 2 ) 4 -> {-1,0}: 



1 ( 66' ^ ( Ul ' U2 ^3,«4) = (1, 1,0, 1), (0, 1,0, 1), | ^ 
(1) 1) 1, 1) 



« /, ,/ x f -1 (W/^O, (Ul, 112,143, «4) = (1,1, ljl)), , n x 

otherwise. 



Then we have 

C Y (3) = j(N 3 + K 3 ), (10) 
A/3 = 1^(2), (11) 

^ 3 = ^pTP^I ^ ^ ^ ^ 

I ^11 ^1 ^2 I B 2 eS 2 Bi€Bi / 2 eQy 2 6e(F 2 )4 

MH B *)fc b )(®i(bi; b) + $ 2 (6i, 6 2 ; 6) + * 3 (&i, fc; 6)), 

where 

0i(/i;6) := tt{ i e [Li] | - + l),/i(i + 2)) = 6}. 

Proof. Equations (I10p and (jilt) have already been shown. From (I17p in Appendix El we find 
that 

1 Ll 

^ 3 = |fiij|S 2 ||fiyJ S Y Y Y 

i ^1 *u Y2 < b 2 £B 2 b^B! f 2 en Y2 i=i 

^ 1 (b i ;I(B 2 )f 2 (i-l),I(B 2 )f 2 (i),I(B 2 )f 2 (i + l),I(B 2 )f 2 (i + 2)) 
+ $2(6,, b l+1 ;I(B 2 )f 2 (i - 1), I(B 2 )f 2 (i), I(B 2 )f 2 (i + l),I(B 2 )f 2 (i + 2)) 
+ $ 3 (6i-i, i"(B2)/2(t - 1), I(B 2 )f 2 (i), I(B 2 )f 2 (i + 1), 7(S 2 )/ 2 (t + 2)) 

Since 

(6i,6 2 , ... ,6^) G Bi — s- (6j,6j + 1, . . . . . . ,6i_i) G Si, 

we can rewrite 

£ $i(bi;I(B 2 )f 2 (i - l),I(B 2 )f 2 (i),I(B 2 )f 2 (i + l),I(B 2 )f 2 (i + 2)) 

= I(B 2 )f 2 (i-l),I(B 2 )f 2 (z),I(B 2 )f 2 (i + l),I(B 2 )f 2 (i + 2)). 

BieBi 

Furthermore 
Li 

^ $i(6r, J(S 2 )/ 2 (i - 1), I(B 2 )f 2 (i),I(B 2 )f 2 (i + 1), J(fl 2 )/ 2 (i + 2)) 

i=l 
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= £ «M/(£ 2 )/ 2 ;b)$i(&i;b). 

6e(F 2 ) 4 

Similar relations hold for <J> 2 an d $3, and we obtain 

Ks = IBillBallfiyJ ^ ^ ^ ^ 

Ml(B 2 )f 2 ; b)($i(6i; b) + $ 2 (6 1; 6 2 ; b) + $3(61, &35 &))• 
Theorem 3. 

Cy(3) = i ^1(7^(2) + 2#J 1) (3k 3 + 2fc 2 - 7^ 2 ) 



+ - 4fc 3 - 3A: 2 - ft 2 - (JT^ + l)L 2 Cy 2 (l) - L x C Yx {h 2, 3)} 

fei - 2 



^)^M^^ j (12) 



where 

r(j -) kj-i + l 



Lj ~t~ /lj — i 
and 

2k 3 k 2 k 2 (k 2 - 1) 



n 2 = - + 



L 2 + fe 2 - 1 L 2 + /c 2 - 2 
Proof. From (|18p and (|19p in Appendix [Bj we obtain the following equations: 

E 2 0i(/(S 2 )/ 2 ;b)$ 1 (6 1 ;b)= [ Ll + ^ 1 1 " 2 ) {2 (3^3 + 2^-^(52))}, (13) 
BieSi be(F 2 ) 4 ^ 1 ' 

E (Ml(B 2 )fc 1, 0, 0, 0)* 2 (6l, 62; 1, 0, 0, 0) + 0i(J(B 2 )/ 2 ; 0, 0, 0, 0)$ 2 (bi, 62; 0, 0, 0, 0) 
S16B1 

+ 4>\{I(B 2 )f 2 \ 1, 1, 0, l)$ 2 (bi, 62; 1, 1,0, 1) + Mi(B 2 )fr, 0, 1, 0, 1)^2(61, 62; 0, 1, 0, 1) 
'L\ + h - 3' 



v >i (i(s 2 )/ 2 ; 1, 0, 0, 0) + </>i{i(B 2 )f 2 ; 0, 0, 0, 0) 

L (I(B 2 )/ 2 ;l,l,O,l)-0i(I(5 2 )/ 2 ;O,l,O,l)) (14) 



= ( Ll , + *\ 3 ) fii - 4fc 3 - 3fc 2 - j(B 2 ) - E C(6i; / 2 (i), / 2 (i + 1)) " E " 3)) , 

^ 1 ' V i=l i=4 / 

E Ml(B 2 )f 2 ; 1, 1, 1, l)($ 2 (bl, 62; 1, 1, 1, 1) + 1, 1, 1, 1)) 

BieBi 



V fei - 2 / V kx-2 
where (b\,b 2 , . . . , 6^ 2 ) = B 2 £ B 2 and 

m \ I 1 (M0, (m,n 2 ) = (l,l)), 

C(6;«i,«2) := < _ . (16) 

U otherwise. 
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Note that, from the definition, 

|g L , E E M^fr, 1,1, 1,1) = 1^(1,2,3). 

1 2,1 y2 ' B 2 eB 2 / 2 eQy 2 

Since 

(bi,b 2 , . . . ,b L2 ) £ B 2 — 5- (6i,6i+i,...,6 i2 ,6i,...,6i_i) € S 2 

and 

H'ir in / 2 G n y j = K'n' m f 2 G OyJ, 

£ 5^C(6.;/a(i),/a(i + i))= E C(&i;U)(x>(i-3)) = ( L2 + h ~ 2 )l 2 c Y2 (i). 

B 2 £B 2 i=l B 2 eB 2 \i=4 / ^ 2 ' 

Hence, we have 

r /Li+fci-2\ 
1 ^ b 2 gb 2 I I fci J 

/ Pi \ /Li+fei-3\ 

+ l Ll _ 4 fc 3 - 3fc 2 - j(5 2 ) - (C(6i; 1, 1) + l) E W ~ 3 ) 

\ i=4 /I fci J 

( L ev 3 ) + (^i-v 4 ) i 



m v \ E 1 (/(^)/ 2 ;i,i,i,i) 

2fci (3fc 3 + 2fc 2 - ft?) 
Li + kx - 1 

fci (fci - 1) {Li - 4fc 3 - 3fc 2 - £ 2 - ( x 2 q^i + l) ^ 2 Cy 2 (l) - LiCV^l, 2,3)} 

+ {Lx+hi - l)(Li + fci -2) 

fci(fci-l)(Li-l)LiC yi (l,2,3) 



(Li + fci - l)(Li + fc x - 2)(Li + fci - 3) 
From Proposition [21 this equation and K^' give (|12p . 



4 Concluding remarks 

In this article, we showed a recursion formula for the iV-point functions Cy(l,2, ... ,N — 1). 
Using the formula, we explicitly calculated the second and the third nearest neighbor correla- 
tion functions. In these estimations, we found that Cy(3) is essentially obtained from Cy(l), 
Cy(2) and Cy(l,2,3). We expect that such recurrence formulae may exist for general iV-point 
functions. Obtaining the recurrence formulae and to clarify their relation to correlation functions 
for quantum integrable systems are problems that will be addressed in the future. 

A Evaluation of 7£ 3 in Cy-(3) 

Let B 2 G £> 2 , f 2 g J7y 2 and 

R(B 2 ;f 2 ) - 7^ E (»{' 1101 '> ' 1011 ' € I(Bi)I(B 2 )f 2 } 
1 11 BieBi 
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- |j{ £ 100r, '111P € hI(B 2 )f 2 disappear in I{B 1 )I(B 2 )f 2 

To evaluate R(B 2 ; f 2 ) we have to count the balance of the increment and the decrement of the 
1 * *1 patterns by 10-insertion. There are two types of such 1 * *1 patterns: 

(A) One of the Is at the ends of a 1 * *1 pattern originally belong to I\I{B 2 )f 2 ] 

(B) Both of the Is at the ends come from 10-insertion I 2 (B\). 

In each case, the variation is listed as follows. 

(A) Let fi = I(B 2 )f 2 . When a block 1010- • ■ 10 is inserted between /i(n) and fi(n + 1), the 



number of 1 * *1 patterns of type (A) is given as: 

00 : 

— ► 001010- •• 1000 => ±0(+0,-0), 

— > i0oh)io---iooo => +l(+l,-0), 

— > 001010- •• 1001 => +l(+l,-0), 

— >• 10OH)1O- • • 1001 =4> +2 (+2,-0); 



(1) /i(n)/i(7i + l) = 

(a) 0000 

(b) 10000 

(c) 0001 

(d) 10001 

(2) /i(n)/i(n + l) = 

(a) 1000 - 

(b) 1001 - 

(3) /i(n)/i(n + l) = 

(a) 0011 

(b) 10011 

(4) / 1 (n)/ 1 (n + l) = 

(a) 01100 

(b) 011100 
(b') 111100 

(c) 0111 

(d) 011111 
(d') 111111 



10 : 
01 



igio io ■ 
lolio 10 ■ 



1000 
1001 



+i (+1,-0), 
+1 (+2,-1); 



001010- •• 1011 
10010 10- •• 1011 



=> +1 (+1, 

=> +1 (+2, 



-o), 
-i); 
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011010- 
0111010 
1111010 
011010- 
0111010 
1111010 



10100 

•10100 

•10100 

1011 

•10111 

•10111 



+2 (+2,-0) 

+1 (+2,-1) 
+1 (+3,-2) 

+1 (+2,-1) 
±0 (+2,-2) 
±0 (+3,-3) 



where +1 (+1, —0) etc. denote the numbers of the total (created, disappeared) 1**1 
patterns. 

Note that total number does not depend on the length of the blocks and that the variation of 
1 * *1 patterns is determined for the four elements f\{n — 2)f\(n — 1) f\{n) f\{n + 1). Hence, for 
b n ^ 0, the contribution from a boundary and a block is summarized as 

' ((0,0,0,0), (1,1,1,1)), 

/ (1,0,0,0), (0,1, 0,0), (1,1, 0,0), (0,0,1,0), 
1 (1,0,1,0), (1,1, 1,0), (0,0,0,1), (0,1, 0,1), 
V (1,1, 0,1), (0,0, 1,1), (1,0, 1,1), (0,1, 1,1) 
k 2 ((0,1,1,0), (1,0,0,1)). 

For example, there is no contribution for (/i(n — 1), f\{n + l),/i(n + 2)) = (0,0,0,0), 
(1,1,1,1), and +2 for (/i(n - l),/i(n),/i(n + l),/i(n + 2)) = (0,1,1,0), (1,0,0,1). 

(B) When blocks are inserted between f\(n — 1) and fi(n) and between f\{n) and f\{n + 1) 
(6 n _i ^ 0, b n ^ 0), the total number of 1 * *1 patterns generated between two blocks is 
given as: 
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(1) / 1 (n-l)/i(n)/i(n + l) = 000: 

— > 010- ■■ 1001010- •• 100 (+1,±0), 

(2) / 1 (n-l)/i(n)/ 1 (n + l) = 001: 

— ► 010- • -10010 10 •••101 (+1,-1), 

(3) / 1 (n-l)/i(n)/i(n + l) = 010: 

— > 010- 101010 10 -100 => (±0,±0), 

(4) / 1 (n-l)/i(n)/i(n + l) = 011: 

— ► 010- ■■ 1011010--- 101 (+1,-1), 

(5) / 1 (n-l)/i(n)/ 1 (n + l) = 100: 

— ► 1010- ■■ 1001010- ■■ 100 => (+1,-1), 

(6) / 1 (n-l)/i(n)/ 1 (n+l) = 101: 

— ► 1010- ■■ 1001010- ■■ 101 => (+1,-2), 

(7) / 1 (n-l)/i(n)/i(n + l) = 110: 

— ► 110---10l|^1010---100 =4> (±0,±0), 

(8) h {n ~ l)/i (n)/i (n + l) = 111: 

(a) — ► 0110 • • • 10110 10 • • • 101 (+1,-1), 

(b) — > 1110- •• 1011010- •• 101 (+1,-2); 

Here (+1,±0) etc. denote the number of 1 * *1 patterns which (arise between blocks, are 
eliminated by the 10-insertion of blocks). 

When a block is inserted between fi(n — 2) and f\(n — 1) and another block is inserted 
between /i(n) and f\{n + 1); 

(9) h{n-l)h(n)h(n + I) = 111: 

(a) — ► 010- ■■10111010--- 101 (±0,±0), 

(b) — ► 110- •• 101110 10 •••101 (±0,-1). 

From these list, we find that the pattern fi(n— l)/i(n)/i(ra+l)/i(n+2) determines the variation 
locally. For b n b n+ i ^ 0, 

1 ((1,0, 0,0), (0,0, 0,0)), 

(1,1, 0,1), (0,1, 0,1), 

(1,1,1,1) 
otherwise, 

and for 6 n _i / 0, 6„ = 0, b n+ i ^ 0, 

f "I (1,1,1,1), 
otherwise. 

Therefore, using the functions $i, $2 and $3 ([Zl+dH]) in Proposition [21 we obtain 

1 11 BieBi t=l 



^ 1 (b i ;I(B 2 )f 2 (i-l),I(B 2 )f 2 (J,),I(B 2 )f 2 (i + l),I(B 2 )f 2 (i + 2)) 
+ $ 2 (6 i5 b i+1 ;I(B 2 )f 2 (i - 1), I(B 2 )f 2 (i),I(B 2 )f 2 (i + 1), /(B 2 )/ 3 (t + 2)) (17) 
+ $ 3 (&i-i, b i+1 ;I(B 2 )f 2 (i - l),I(B 2 )f 2 (i),I(B 2 )f 2 (i + 1), 7(S 2 )/ 2 (i + 2))) . 
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B Derivation of equations (1131). (1141) and (1151) 

From (|17p and Proposition we find 

R{B 2 ;f 2 ) = ^- E Ml(B2)f2;b)(^ 1 (b 1 ;b) + ^ 2 (b 1 ,b 2 ;b) + ^ 3 (b 1 ,b 3 ;b)). 

1 11 BiSBi 6G(F 2 ) 4 

Let fi := I(B 2 )f 2 . First we consider the term: 

Ra(B 2 -J 2 ):=-±- E E MlWfaQQtfcb) 

' 1 ' BieBi 6e(F 2 ) 4 

= lbi E E $ i( &i '"/ i ( n - 1 )'/ i ( n )'/ i ( n + 1 )'/ i ( ri + 2 ))- 

' X ' BieBi n=l 

For 6 / 0, we observe $i(&;b) = #{'10701' € b} - #{'101', '010' € b}. There are k 2 '10's, the 
same number of '01's and j{B 2 ) + j|{'10' inserted between consecutive l's} '101's in f x . Since 
the number of '010's is equal to jj{'10'} — ft{'10' inserted between consecutive l's}, we find 

Ra(B 2 ; f 2 ) = 6k 2 - 2j(B 2 ) - 2k 2 . (18) 
The rest terms can be evaluated by counting the numbers of 3-tuples in f x . 

• the number of (f x (n - 1), f x (n), f x (n + 1)) = (0,0,1), (1,0,0) is 

2fc 2 -jt{n€[L 2 ]|(/i(n-l),/i(n),/ 1 (n + l)) = (1,0,1)}; 

• the number of (/i(n — 1), fi(n), f x (n + 1)) = (0, 1,0) is 

L 2 

i=l 

where ((b;u x ,u 2 ) is defined by (fl~6|) : 

• the number of (/i(n - 1), /i(n), f x (n + 1)) = (0, 1, 1), (1, 1, 0) is 

L 2 

2fc 3 + 2^c&;/2«,/2(*' + i)); 

i=l 

• the number of (f x (n - 1), /i(n), f x (n + 1)) = (1, 0, 1), (1, 1, 1) is 

Pi 

J2h(i-3)+j(B 2 ). 

Hence we have 

tt {'000' € I(B 2 )f 2 } - (t {'101' G I(B 2 )f 2 } 

L 2 Pi 

= L X - ik 3 - 3k 2 - j(B 2 ) - CiW, / 3 (»), f2(i + 1)) - E - 3)- 

j=l i=4 

Therefore 

cf> x (I(B 2 )f 2 ; 1, 0, 0, 0) + MHB 2 )f 2 ; 0, 0, 0, 0) - <t> x {I{B 2 )f 2 ; 1, 1,0, 1) 

- 0i(iOB 2 )/ 2 ; o, l, o, l) = tt {'ooo' g /(£? 2 )/ 2 } - tt {'ior g i(s 2 )/ 2 } 

L 2 Pi 

= L x - ik 3 - 3k 2 - j(B 2 ) - CiW, h(i), f2(i + 1)) - E " 3 )- ( 19 ) 

i=l i=4 

From (fl8j) and (HH), R(B 2 ; f 2 ) can be obtained as shown in the proof of Theorem [3l 



16 



J. Mada and T. Tokihiro 



Acknowledgements 

The authors wish to thank Professor Atsushi Nagai for useful comments. 

References 

Takahashi D., Satsuma J., A soliton cellular automaton, J. Phya. Soc. Japan 59 (1990), 3514-3519. 

Yura F., Tokihiro T., On a periodic soliton cellular automaton, \J. Phys. A: Math. Gen. | 35 (2002), 3787- 
3801, |nlin.SI/011204l) 



Fukuda K., Okado M., Yamada Y., Energy functions in box ball systems, Internat. J. Modern Phys.~%\ 15 
(2000), 1379-1392, |math.QA/9908116l 

Hatayama G., Hikami K., Inoue R., Kuniba A., Takagi T., Tokihiro T., The A^f automata related to 
crystals of symmetric tensors, \J~Math. Phys] 42 (2001), 274-308, |math.QA/9912209j 

Mada J., Tokihiro T., Correlation functions for a periodic box-ball system, | J. Phys. A: M ath. Theor\ 43 
(2010), 135205, 15 pages. larXiv:0911.3953l 

Jimbo M., Miwa T., Quantum KZ equation with |q| = 1 and correlation functions of the XXZ model in 
the gapless regime, \J7Phys. A: Math. Gen. \ 29 (1996), 2923-2958. 

Mada J., Idzumi M., Tokihiro T., The exact correspondence between conserved quantities of a periodic box- 
ball system and string solutions of the Bethe ansatz equations, \ J. Math. Phys. 47 (2006), 053507, 18 pages. 

Torii M., Takahashi D., Satsuma J., Combinatorial representation of invariants of a soliton cellular automa- 
ton, \Phyr]J\ 92 (1990), 209-220. 



